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Nature versus Nurture: The curved spine of the galaxy 
cluster X-ray luminosity — temperature relation 
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ABSTRACT 

The physical processes that define the spine of the galaxy cluster X-ray luminosity 
- temperature (L-T) relation are investigated using a large hydrodynamical simulation 
of the Universe. This simulation models the same volume and phases as the Millen- 
nium Simulation and has a linear extent of 500/i~^Mpc. We demonstrate that mergers 
typically boost a cluster along but also slightly below the L-T relation. Due to this 
boost we expect that all of the very brightest clusters will be near the peak of a merger. 
Objects from near the top of the L-T relation tend to have assembled much of their 
mass earlier than an average halo of similar final mass. Conversely, objects from the 
bottom of the relation are often experiencing an ongoing or recent merger. 
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1 INTRODUCTION 

Since the launch of the XMM-Newton and Cha ndra satel- 
lites l|jansen et al.l I2001I : IWeisskopf eraLll2002l ). measure- 
ments of the X-ray emission from hot gas in clusters of 
galaxies have achieved unprecedented levels of accuracy and 
depth. However, the physical origin of the scaling relations 
between observable quantities, such as the luminosity of the 
X-ray emitting gas and its temperature, remain only partly 
understood. 

T here are currently a n umber of surveys (iRomer et al.l 

I2OOII : ISchwope et al. I I2OO3I : iPierre et al. I I2OO6I ) in progress 
with the potential to greatly expand our understanding of 
the processes that define correlations such as the luminosity- 
temperature (L-T) relation of clusters. For this potential to 
be realised we require a sound theoretical basis upon which 
to work. To this end, numerical hydrodynamical simula- 
tions have become indispensable tools and continue to grow 
in size and complex ity (Pearcc ct al. 2000; Kav ct al. 2007; 
iFaltenbacher et al][2007l ) but they have to date lacked a suf- 
ficiently large dynamic range in mass. In this work we use a 
hydrodynamical model of a large volume that contains over 
a hundred galaxy clusters. For the first time we are able to 
study the evolutionary processes within a cosmological con- 
text as we have hundreds of well resolved objects spanning 
a large dynamic range rather th an the more typical handful 
jRowlev Thomas, fc Kavll2004l) (hereafter R04 ), or idealised 
models ^Ritchie fc Thomasll2002l : iPoole et al.ll2006l ) . 



This paper is organised as follows; in the remainder of 
this section we summarise the work done to date on defin- 
ing the physical processes that define the shape of the L-T 
relation. Then, in section 2, we give an account of the sim- 
ulations we have undertaken, explain how our cluster sub- 
sample was selected and how the properties of these clusters 
were derived. Section 3 details our results before we discuss 
their implications and conclude in section 4. 

X-rays are chiefly emitted from the hot gas in clus- 
ters via thermal bremsstrahlung (for dark matter halos 
more massive than lO^^/i"^]^© their temperature is typi- 
cally above 2keV). For such a homologous population lKaiserl 
(1986) showed that simple scaling relations were expected 
between bulk properties such as the mass, temperature and 
luminosity. Observational work subsequently found that the 
properties of X-ray clusters where indeed related but the 
slopes of the relations were not those derived by Kaiser. 

iKaised (|l986l ) assumed that galaxy clusters were self- 
similar entities and that therefore only a single property, 
such as the mass, was required in order to describe the 
other bulk properties. Such a homology results in an L- 
T relation with a power-law slope of 2. However, as fig- 
ure [1] demonstrates, X-ray observations of clusters with a 
median redshift of ~ 0.07 found that the slope was closer 
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Figure 1. Compilation of low-redshift observed group and cluster X-ray luminosities within rsoo compared to their 
emission-weighted temperature, rsoo is a radius enclosing an overdensity of 5 00. The sma ll points arc the simulated 
group s and clusters used in this work. The data was taken from var i ously: | Markev it"ch^ (1998); Arnaud & EvrardI 
('l999''):'Wu et al.' ('l999''):l HeIsdon fc PonmanI l|2000t l: IXue fc Wul (|200ol '): lHorner et al.l l|2001l '): IMuIchaev et aLipOof h 
FOsmond fc Ponman l20ol) ! 



I2OO3I : lOsmond fc Ponmari[ 2004 ) and perhaps bec ame even 
steeper on group scales ( Helsdon fc Ponmanll2000l l. 



Hydrodynamical simulations performed in the absence 
of cooling or any additional heat sources other than com- 
pression and shock heating have lon g been known to repro- 
duce the self-similar hierarchy well ([Eke. Navarro, fc FrenkI 
1 19981 : iNavarro. Frenk. fc White! 1 19951 ). Unfortunately they 
do not reproduce either the slope or the normalisation of 
the observations, producing clusters that are too bright for 
any given temperature, even at the bright end. Following 
this work, simulations with limited physics within a cos- 
mological volume have been used in an attempt to recon- 
cile the apparent discrepancy between theory and observa- 
tion regarding the slop e of the L-T relation ([P earcc et al.l 
2000; MuanwouE ct al. 200ll: iBialek. Evrard. fc Mohill200ll ; 
Borgani et al. 2002 ) . Those models showed that a simple 



cooling or preheating scheme was sufficient to match the 
simulated L-T relation t o that observed at redshift zero. 
More recently iKav et al] (|2007l) investigated the effects of 
feedback on the X-ray properties of clusters in hydrodynam- 
ical simulations, and demonstrated that their results were 
in good agreement with both the observed scaling relations 
and structural properties (e.g. entropy and temperature pro- 
files), particularly for cool-core clusters. 

iBalogh et al.l (|2006l ) investigated the role that preheat- 
ing, cooling and concentration of the halo profile can have on 



the scaling relations. They found that, for a realistic range 
of halo concentrations, the scatter generated was minimal in 
comparison with observed values. Variations in the cooling 
time of the gas in the centre of clusters could account for 
much of the scatter but is limited by the age of the universe 
and so could not explain the whole range. Finally, varying 
feedback from supernovae and AGN could explain the en- 
tire range, but required an order of magnitude difference in 
energy injection to cover the whole envelope. Their result 
implies that it is processes in the cores of clusters that are 
primarily responsible for driving the sc atter in the scaling re - 
lati ons. Thi s confi rms e arlier work by iFabian et al. (Il994t ): 
.MarkevitchI (| 19981 ') and iMcCarthv et al. (2004). K ay et al I 
(|2007^ identify the scatter with strong cool core clusters, 
and expect the scatter to be smaller at high redshift due 
to the diminished prevalence of such systems. Nowadays, 
the general consensus is that the scatter is largely due to 
the strength of the X-ray core. In this work, which includes 
strong preheating, X-ray cores are absent. This allows us 
to study the shape of the relation without the additional 
complication of a large intrinsic scatter. 

In this work we will use a sample of halos identified from 
the full simulation volume. With these we will show that be- 
cause mergers tend to move clusters up the L-T relation they 
extend it beyond the point where the most massive, relaxed 
clusters are expected to lie. Thus many of the brightest, most 
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luminous objects are ongoing or recent merger events which 
(as R04 point out) may be difficult to resolve observation- 
ally if they are close to the peak of the merger. In addition, 
because we have many closely spaced outputs we can track 
the motion of each of our clusters on the L-T plane, allowing 
us to define a "mean merger" vector. As this vector is not 
perfectly parallel to the L-T relation but rather falls slightly 
below it, a gentle roll in the relation naturally arises. 



2 THE SIMULATIONS 

The simulation used i n this work i s part of the Millen- 
nium Gas Simulations (|Pearce et al.1 ,2007'). This sequence 
of hydrodynamical simulation s all have the same volume 
as the Millennium Simulation l|Springel et al.l 120051 ) as well 
as utilising the same amplitude and phase for the ini- 
tial perturbations. The cosmological parameters for both 
the Millennium Simulation and the gas counterparts were: 
Oa = 0.75, f^M = 0.25,^6 = 0.045, /i = 0.73, n = 1, and 
as = 0.9, where the Hubble constant is characterised as 
100/ikms~^Mpc~^. These cosmological parameters are con- 
sistent with recent combined analyses from WMAP data 
Spergel et ail 120031) and the 2dF galaxy redshift survey 



CoUess et al.ll200l[ ). The simulation volume is a comoving 



cube of linear size 500/i ^Mpc containing 500 million dark 
matter particles and 500 million gas particles. Their masses 
are 1.422 x lO^°/i"^M0 and 3.12 x l(fh-^MQ respectively. 
The simulation includes radiative cooling of the gas, with the 
metalicity set at a constant value of O.S Zq, similar to that 
observed within the intra-cluster medium (jSarazi nll986D and 
pr eheating. The prehea ting is implemented in a similar way 
to lBorgani et al.l (|2004l ): at redshift 4 the whole volume is 
heated to 200keV/cm^, a value chosen such that the result- 
ing L-T relation at redshift matches observations. Star 
particles are formed from cold, dense gas particles when 
a temperature threshold (2 x 10^ K), a density threshold 
(uh = 4.18 X 10~^^ g.cm""^ and an overdensity threshold 
(100 times critical) are all passed, but the process of con- 
verting gas to stars has no effect on the thermal dynamics 
of the system, other than to make the particles coUisionless. 
The effect of the preheating in this simulation is so extreme 
as to prevent any further star formation since redshift 4. 



2.1 Sample selection 

At redshift zero the entire volume was processed to ob- 
tain a set of friends-of-friends halos with a linking length 
of 0.2 times the mean interparticle separation. Within each 
of these halos the most bound particle was found and used 
as the centre for a spherical over-density calculation that 
extended to r2oo, a radius enclosing an over-density of 200 
times the cosmic mean. The analysis presented in this work 
is for a fixed radius of rsoo (the radius enclosing an over- 
density 500 times the the cosmic mean density), roughly 
0.59 X raoo for NFW halos l|Navarro. Frenk. fc Wliit3ll9971 ) 
of typical concentration. Within this radius we calculate the 
bolometric luminosity an d the emission weighted tem pera- 
ture assuming a standard (ISutherland fc Dopita|[l993l ) cool- 
ing function for a uniform metalicity gas of O.3Z0. 

Three sub-samples were selected from the top, bottom 
and median of the L-T relation. We refer to the sample of 
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Figure 2. The clusters selected for further analysis: the high 
and low scatter samples are shown as open circles and squares 
respectively, above and below the line which indicates our fitted 
moan relation. The control clusters are shown as solid circles. The 
underlying faint points show the full sample at redshift 0. 



clusters that are more luminous and/or cooler than expected 
as coming from the top of the relation, with conversely under 
luminous, hot clusters coming from the bottom. We also 
select a control sample of clusters from close to the median of 
the relation. The clusters were selected such that the range 
of masses within each sample spanned the entirety of the 
available relation. We only consider objects containing more 
than 1000 particles and that are at least two virial radii 
away from any larger neighbour. This ensures a meaningful 
estimate of the cluster bulk properties. 

Once selected at redshift zero, each of our 108 clusters 
was traced backwards in time until their mass dropped below 
our imposed resolution threshold of 1000 particles. The final 
locations of the selected clusters on the L-T plane are shown 
on figure (2] where the high scatter clusters are indicated by 
open circles, low scatter by open squares and the control 
sample by filled circles. The full sample is shown faintly in 
the background, together with our fitted median relation 
indicated by the line. Once the mass accretion histories of 
the sample had been extracted a small amount of smoothing 
was introduced in order to remove merger induced ringing 
in the cluster mass. 



3 HALO PROPERTIES 
3.1 Low-scattered halos 

The L-T histories of nine low scattered halos ranked by final 
mass are shown in figure O with their respective mass ac- 
cretion histories given in figure |4l The dark track in the L-T 
plane follows the history of each object from redshift 1.5 to 
the present day. The background points show the location 
of the entire sample of clusters in the L-T plane at z = 0. To 
the bottom right of each panel are two additional vectors. 
The short line shows the mean evolution of the control sam- 
ple between z = 0.5 and the present. The other longer line 
shows the evolution of each particular group over the same 
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Figure 3. X-ray evolution in the L-T plane of 9 of the low-scattered clusters ranked by final mass. In each of the nine frames the L-T 
relation at redshift zero is defined by a sample of clusters (background points) whilst the track showing the evolution of each cluster is 
shown as the dark line. To the bottom right of each panel the evolution in the L-T plane for each cluster (long dashed line) and the 
mean of the control sample members of similar mass (short solid line) between a redshift of 0.5 and 0. 



redshift interval. While the control sample moves slightly up 
in luminosity and down in temperature, 8 of the 9 low scat- 
tered clusters move dramatically to larger luminosities and 
temperatures, nearly parallel to the spine of the L-T rela- 
tion. This is in agreement with the trend noticed by R04. 
Figure |4l which shows the corresponding mass accretion his- 
tories for these objects demonstrates that 8 of the 9 low scat- 
tered objects are in the process of an ongoing major merger 
and are much hotter and brighter than would be typical for 
objects of their mass. In each of these panels merger events 
are denoted by the bold sections of line. The masses at red- 
shift 1.5 are significantly lower than expected in all bar one 
case, (the mean mass accretion history for objects of each 
mass is indicated by the dotted line on figure |4)| . 

3.2 High-scattered halos 

Figures [5] & |6] show the L-T evolution and mass accretion 
histories for nine of the high scattered clusters. These nine 
halos end up significantly above the mean relation and as 
can be seen from their L-T tracks in figure [5] eight of the 
nine (all except panel e) slightly lose temperature rather 
than gain temperature along with the mean of the control 
sample. The mass accretion history makes it clear why this 
is the case: all bar panel (e) assemble their final mass early, 
with significantly more mass in place at z = 1.5 than that 
collected by the control sample. The object in panel (e) has 



just undergone a merger. We conclude, as did'Balog h et al] 
(2006), that high-scattered objects are in general early form- 
ing with consequently slightly more concentrated dark mat- 
ter profiles resulting in slightly more luminous objects at a 
particular final mass. 

3.3 Properties of mergers 

As discussed in the previous section the motion of an ob- 
ject on the L-T plane during a merger is a significant driver 
behind finding it below the mean L-T relation at any given 
mass, particular at the high-mass end. To explore this fur- 
ther we extracted a sample of mergers from the mass accre- 
tion histories of the clusters used previously in this work. 
In order to distinguish a merger from gradual accretion we 
require that a cluster gains significant extra mass over a 
short period of time. Specifically we define a merger in the 
following way: 

• A growth in mass through the merger event of at least 
15% of the cluster's final mass. 

• A ratio of at least 1:4/3 between the mass before the 
merger and the mass at the peak of the merger. 

• The mass accretion rate must exceed 14% of the final 
mass per Gyr. 

Merger events automatically identified using this proce- 
dure are shown on the mass accretion history figures |4] & [6] 
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Figure 4. Mass accretion tiistory for tiie 9 low-scattered clusters whose L-T evolution was shown in figure[3] The lighter line is the mass 
accretion history for the mean of similar mass clusters in the control sample whereas the darker line shows the history of the individual 
cluster. Bold sections demote periods when the cluster is undergoing a merger as defined in section 3.3 below. The mass of the cluster 
at each expansion factor is normalised by its final mass, labelled as "scaled mass" . Also plotted is a horizontal dotted line to show when 
a cluster has assembled 70% of its final mass. 



as the bold sections of the lines. The peak of a merger is con- 
sidered to be the point at which the cluster's mass is great- 
est. As the dark matter halos subsequently pass through 
each other the final mass is usually below this value. As can 
immediately be guessed by simply comparing the number 
of bold line sections in figures |3] & [5] the mean number of 
mergers undergone by clusters in the low-scattered sample 
is over three times higher than clusters in the high-scattered 
sample. 

By identifying the location on the L-T plane of each 
object at the start and peak of each merger we can produce 
a "cricket score" diagram (figure [7]), where each line denotes 
the motion on the L-T plane due to one merger. As the 
merger timescale is short the net drift of the relation is small 
while the merger is ongoing. As can be clearly seen, the net 
effect of a typical merger is to move an object up the relation 
and on average slightly below it. This tendency for mergers 
to fall below the mean relation is further evidenced by the 
large number of ongoing mergers present amongst the low- 
scattered objects. 

Interestingly, an "average merger" vector, indicated by 
the dotted line on figure [T] closely parallels the slope of the 
very high mass clusters. The tendency for mergers to boost 
clusters along, but at a slight angle to, the relation also 
drives a slight roll that is found at the high mass end of the 
L-T relation. 



4 DISCUSSION 

This work examines the physics that underlies the spine 
of the X-ray L-T relation. Due to our strong preheating 
prescription our halos do not have strong cores and as 
such do not reproduce the large scatter in the observed 
L-T relation, allowing us a clear window into the basic 
physics. We intend to examine the physical origin of the 
observed scatter in future work (iGazzola et akl uOOT^) where 
a more physically motivated energy feedback prescription 
will be used and bright cooling cores are present. Preheat- 
ing schemes such as the one used here are well known to 
accurately reproduce t he slope and normalisation of the L- 
T relation as a whole p^earce et al.l 20001: iMuanwong et al.l 



I2OOII : iBialek. Evrard. fc Mohd l200ll : iBorgani et al.ri2002l '). 
The model we have implemented also accurately reproduces 
the mean location of halos on the L-T plane at the present 
day but in a much larger volume than has typically been 
used previously. In the real world bright cooling cores will 
further complicate matters but the processes discussed here 
which relate to the outer halo properties will underlie these, 
with the variation in core properties leading to a scatter 
about the relation discussed here. 

By identifying mergers using the mass accretion histo- 
ries of our objects and matching these episodes to the motion 
of each object on the L-T plane we have derived a "mean 
merger" vector in this plane. This vector lies largely parallel 
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Figure 5. L-T evolution for 9 of the high-scattered sample. The symbols and lines have the same meaning as in figure |3] 



to the cluster L-T relation, as previously noted by R04. At 
any particular time the mass function of the dark matter 
haloes present within a volume will be exponentially trun- 
cated at the high mass end above some characteristic mass 
scale. The large boost generated during a merger will pro- 
duce points on the L-T plane appearing to lie above this 
characteristic mass, where there should be few objects. We 
therefore expect the majority of the brightest objects to be 
experiencing ongoing mergers, although they may be diffi- 
cult to identify if they are close to their peak. 

The mean merger vector we have derived is not exactly 
parallel to the L-T relation but rather lies slightly below 
it. This behaviour leads to all bar one of our low-scattered 
objects being obvious recent or ongoing merger events (fig- 
ure |4]). We also note that at the high mass end the vast 
majority of our haloes lie below the mean relation shown on 
figure [2] The fact that the mean merger vector lies slightly 
below the mean relation provides a natural explanation for 
the slight curvature evidenced in the simulated relation. 

In summary, while it is straightforward to reproduce the 
observed slope and normalisation of the X-ray luminosity- 
temperature relation using a simple preheating scheme, such 
a scheme does not reproduce the observed scatter. As a pre- 
heating model includes the full underlying framework of the 
hierarchical build up of structure bulk mergers are not signif- 
icant drivers of this scatter. Mergers can, however, produce 
objects that are brighter and hotter than would be expected 
from the cluster mass as merger events drive objects along 
the L-T relation towards the bright end. We find that a typ- 
ical merger track does not exactly parallel the L-T relation 
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Figure 7. Relative motion on the L— T plane during each of the 
mergers defined in section 3.3. Each line represents a single merg- 
ing event. The long dashed line indicates the mean L— T relation 
whereas the dotted line indicates the mean merger direction. 



but rather lies slightly below it, leading to a prevalence of 
recent or ongoing merger events on the low-scattered side of 
the relation. This process also leads to a slight curvature of 
the mean relation at the high mass end. 
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